PILOT-AIDED CHANNEL ESTIMATION FOR OFDM IN WIRELESS SYSTEMS 

5 

RELATED APPLICATIONS 

The present regular utility patent application claims priority to U.S. provisional patent 
application no. 60/1 14,571, filed December 31, 1998. 

10 TECHNICAL FIELD 

The present invention is directed to a method and apparatus for pilot-aided channel 
estimation for orthogonal frequency division multiplexed (OFDM) wireless communication 
systems. In particular, it concerns pilot-aided channel estimation techniques employing 
Fourier transforms and filtering at the receiver. 

15 

BACKGROUND OF THE INVENTION 

Orthogonal frequency division multiplexing (OFDM) is a potential technique to 
transmit high-bit-rate data over indoor and outdoor wireless communication systems. 

20 Channel parameter estimation is an important technique to further improve data transmission 
performance. There are two different types of channel parameter estimators: decision- 
directed and pilot-symbol-aided. 

Decision-directed parameter estimators have been used in OFDM systems for 
coherent detection, co-channel interference suppression and transmitter diversity. For 

25 systems that are noise limited or with little co-channel interference, coherent demodulation 
with estimated channel parameters can be used instead of differential demodulation, and can 
achieve a 2.5-3.0 dB signal-to-noise ratio (SNR) gain. Moreover, for systems with receiver 
diversity, maximal ratio diversity combining (MR-DC), which is equivalent to minimum 
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Mean-square error diversity combining(MMSE-DC) in the absence of co-channel 
interference, can be achieved directly using the estimated channel parameters. For systems 
with co-channel interference, the coefficients for the MMSE-DC can be calculated from the 
estimated channel parameters and the instantaneous correlation of the signals from each 
5 receiver. For OFDM systems with transmitter diversity using space-time codes, channel 

estimation can provide the necessary information required for the decoding of the space-time 
codes. 

Pilot-symbol-aided parameter estimation is another type of estimation approach. Pilot- 
symbol-aided channel estimation for single carrier systems with fading channels has been 

10 analyzed. Both 1 -dimensional (1-D), double 1-D, and 2-dimensional (2-D) filtering 

algorithms have been proposed for pilot-symbol-aided parameter estimation for OFDM 
systems in terrestrial audio and television broadcasting, fixed and mobile wireless 
communications. However, the filtering algorithms for pilot-symbol-aided parameter 
estimation require channel statistics, such as, the delay profile and the Doppler frequency, 

15 which are usually unknown in wireless environments. Prior art patents have suggested time 
or frequency domain correlation using pilot symbols. 

U.S. Patent No. 5,692,015 to Higashi et al. is directed to pseudo-pilot-symbol aided 
channel estimation in a single carrier system. At the receiver, the channel transfer function 
associated with each pseudo-pilot symbol is estimated and a time-domain interpolation is 

20 carried out to determine the remaining channel transfer functions. 

U.S. Patent No. 5,809,083 to Wright discloses the use of a periodically transmitted 
pair of pilot symbols spaced apart by a fixed difference in a single-carrier system. The pilot 
symbols are used to synchronize the transmitter and receiver and are also used in channel 
estimation. Only the pair of pilot symbols is processed by the channel estimator. 

25 Specifically, the channel estimator compares the current channel-impaired received pilot 
symbols with their expected amplitude and phase, and the average amplitude and phase 
distortion is used to estimate the effects of the channel. Average estimates obtained over 
multiple such periods are then interpolated or otherwise filtered to generate symbol-specific 
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amplitude and phase adjustments to apply to the channel-impaired data symbols. However, 
this technique does not address the issue of OFDM communication and does not account for 
the simultaneous effect of Doppler and multipathing. 

U.S. Patent No. 5,307,376 to Castelain et al. is directed to pilot-aided channel 
5 parameter estimation in an OFDM system. Pilot symbols are broadcast using higher power 
than the remaining symbols in a data stream. At the receiver, the noisy channel parameters 
are subject to an inverse DFT, weighting and thresholding so as to focus on the pilot symbols, 
and the result is then subjected to a DFT. While Castelain '376 provides frequency domain 
interpolation, it does not perform this in the time-domain. Therefore, this technique does not 
10 provide channel estimates which are robust to both Doppler and delay profiles due to 
multipathing. 

SUMMARY OF THE INVENTION 

The present invention is directed to realizing robust parameter estimation for a 

15 channel in a synchronized OFDM communication system which is robust to both Doppler 
and multipathing. At the transmitter, pilot symbols are selectively inserted in predetermined 
frequency bins (tones) of each of a number of OFDM blocks, and a plurality of such OFDM 
blocks are sequentially broadcast. At a receiver antenna, the plurality of OFDM blocks are 
processed, and the resulting signal used to form a noisy two-dimensional channel estimate, 

20 one dimension corresponding to a frequency index of the OFDM tones within each block and 
a second dimension corresponding to a time index of each OFDM block. The noisy two- 
dimensional channel estimate is then subject to two 2-D Fourier transforms, and a diamond- 
shaped filter between them, thereby producing a noise-reduced channel estimate. 

An additional facet of the present invention is that it may be adapted for use in mobile 

25 or diversity receivers having a plurality of receiving antennas. In such case, signals received 
at the various antennas are used to form a single channel estimate which is then used to 
combine the signals to arrive at an optimum estimate of the transmitted signal. 
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BRIEF DESCRIPTION OF THE DRAWINGS . 

The present invention can better be understood through the attached figures in which: 
Fig. 1 is a block diagram of a transmitter in accordance with the present invention; 
Fig. 2 is a block diagram of a receiver in accordance with the present invention; 
5 Fig. 3 is a block diagram of a channel estimator in accordance with the present 

invention; 

Fig. 4a shows a sample pilot symbol grid used for a system in accordance with the 
present invention; 

Fig. 4b shows the contour of a 2-D spectrum for the sample pilot symbol grid of Fig. 

10 4a; 

Figs. 5a and 5b compare the effects on word error rate (WER) and the mean squared 
error (MSE) of different reference signals. 

Figs. 6a and 6b compare the effects on the WER and the MSE of using initial and the 
enhanced estimates for different channel models under a 40 Hz Doppler condition. 
1 5 Figs. 7a and 7b compare the effects on the WER and the MSE of using initial and the 

enhanced estimates for different channel models under a 200 Hz Doppler condition. 

Figs. 8a and 8b shows the effects on the WER of using different Doppler shifts for a 
TU channel model (Fig. 8a) and a HT channel model (Fig. 8b). 

20 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

The aforementioned U.S. patent nos. 5,307,376, 5 5 692 5 015 and 5,809,083 are 
incorporated by reference herein. 

Fig. 1 shows an OFDM transmitter 100 in accordance with the present invention. A 
signal b[n,k] in which the indices k and n represent the k lh tone of the n lh block of an OFDM 
25 sequence, respectively, is presented to an encoder 102. The encoded signal a[n,k] is then 
presented to the modulator 104. In the preferred embodiment, since the phase of each tone 
can be obtained by the estimator, coherent phase-shift keying (PSK) modulation is used to 
enhance the system performance. The modulated signal x[n,k] is then subject to an inverse 
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fast Fourier transform (IFFT) 106 before being broadcast by an antenna 108, in accordance 
with well-established techniques for OFDM transmission protocol. 

Before it is broadcast, signal x[n,k] has embedded therein a plurality of pilot symbols 
which have been added by a pilot symbol insertion block 110. Pilot symbol insertion may 
5 take place just before modulation 1 04, as indicated by the solid line 1 10a, or it may take place 
just before encoding 102, as indicated by the broken line 1 10b. Regardless of the particular 
stage within the transmitter that pilot symbols are inserted, the positions of the pilot symbols 
in the resulting sequence of OFDM blocks is important in the present invention, as discussed 
further below. 

10 Fig. 2 shows a receiver 120 in accordance with the present invention. The receiver 

120 is provided with a plurality of antennas 122a, 122b. The received signal at antennas 122a, 
122b is subjected to a fast Fourier transform (FFT) 124a, 124b to thereby produce respective 
channel-corrupted signals xi[n,k], x 2 [n,k] on lines 125a, 125b. While two antennas and two 
signals are shown, it should be kept in mind that any number p of such antennas may be 

15 present in a mobile diversity receiver. It should also be noted that a channel estimator in 
accordance with the present invention will work even if only a single antenna is used. 

The received signals x m [n,k] are sent on to a channel estimator 126 to estimate 
channel parameters. The estimated channel parameters, in turn, are used in combination with 
the received signals in a diversity combiner 128 to form a combined signal y[n,k]. Combined 

20 signal y[n,k] is then subjected to a thresholder 130 to thereby produce an estimate of the 
demodulated signal a[n,k]. The estimate of the demodulated signal a[n } k] is input to a 
decoder 1 32 to thereby produce an estimate of the decoded signal b[n, k] .The estimate of 
the demodulated signal a[n,k] is also input to the channel estimator to refine the estimate of 

the combined signal y[n,k]. Alternatively, in the event that the pilot symbols were inserted 
25 before encoding at the transmitter 100, the decoded signal b[n,k] may be used in place of the 

estimate of the demodulated signal, as indicated by the broken line in Fig. 2. 
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In general, for a diversity receiver such as that shown in Fig. 2, the signal from the 
antenna at the k th tone of the n* block can be expressed as: 

xJn*] = H m [n,k)a[n,k] + w m [n,k] (Eq. 1) 

5 In the above expression, w m [n,k] is additive Gaussian noise from the mth antenna (m = 1, 2, 
3, . . p) at the k th tone of the n lh OFDM block having zero-mean and variance p. It is also 
assumed that w m [n,k]=s are independent for different n=s, k=s and m=s, i.e., for each 
antenna, the noise has the same statistics at all tones k in each OFDM block n. H m [n,k] is the 
frequency response at the k th tone of the n th block corresponding to the m th antenna, and is 
10 assumed to be independent for different m=s but with the same statistics. Finally, a[n,k] is 
the signal modulating the k th tone during the n* block and is assumed to have unit variance 
and be independent for different k=s and n=s. 

If the transmitted symbols for some k=s and n=s are known at the receiver, the 
channel parameters corresponding to these k=s and n=s ("initial channel parameters") can be 
15 estimated by: 

H m [n, k]A x m [n, k)a [n, k] = H[n, k] + w m [n, k] (Eq 2 j 

where w m [n,k]~ w m [n, k] a* [n, k] . 

From the temporal estimation of the channel parameters for some k=s and n=s, a 
20 robust interpolator can be derived to estimate the channel parameters for all k=s and n=s, 
even if the channel statistics are unknown. 

The receiver 120 further includes a channel estimator 140 which first calculates the 
noisy channel estimates H m [n,k] and further processes the noisy channel estimates to arrive 
at a noise-reduced channel estimate H m [n,k] in accordance with the present invention. 
25 Fig. 3 shows a block diagram of a channel estimator 140 in accordance with the 

present invention. The channel estimator includes a multiplier 142 which calculates the noisy 
channel estimate H m [n,k] from a received signal x m [n,k] on line 142a and a selected 

reference signal a*[n,k] received on line 142b from switch 144. When the receiver is 
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processing the first sequence of OFDM blocks, the reference signal is selected to be the 
known pilot symbols for certain n=s and k=s transmitted by the transmitter in this 
synchronized communication system. In such case, the known pilot symbols are retrieved 
from a local pilot symbol memory 146 associated with the receiver. When the receiver is 
5 processing a subsequent sequence of OFDM blocks, the reference signal is the estimate of the 
demodulated signal a[n,k] or, alternatively, the estimate of the decoded signal b[n,k] . It 

should be understood here that multiplier 142 and switch 144 are preferably implemented in 
software. 

The noisy channel estimate H m [n,k] is sent on to a delay and formatting device 148. 

10 Here, noisy channel estimates for all K frequencies of N successive OFDM blocks are 

assembled into a two-dimensional array G m [n l} n 2 ] wherein r\\ denotes the frequency index 

and r\2 the time index. In the array, each OFDM block occupies a single row with the 
corresponding frequencies being aligned along the columns with the row number of the 
OFDM blocks corresponding to the time sequence in which the OFDM blocks were received. 

15 The formatter 148 then zero-pads both the rows and the columns, if needed, so that each has a 
number of elements that is a power of two. 

The two-dimensional array G m [n lt n 2 ] is then sent on to a 2-D IFFT unit 1 50 where a 
two-dimensional transform is performed to form G m [f lt f 2 ], The output of the IFFT unit 150 
is sent on for frequency domain filtering by multiplying by a filter 1 52 to form G m [/j ,f 2 ] . 

20 The output of the filter 152 is then sent on to a 2-D FFT unit 154 to form G m [n lf n 2 ] from 
which one can extract the noise-reduced channel estimates H m [n, k] . 

Since the 2-D spectrum is usually of a diamond shape, a diamond shape filter is 
preferably used to filter out the additive noise which often occupies the 2-D region 
corresponding to high Doppler and high time delay. The area of the diamond is designed to 
25 match the maximum delay spread and Doppler frequency. Consequently, the estimation 
performance is insensitive to different delay profiles and Doppler frequencies, and so 
provides a robust estimator. 
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In view of the foregoing, it is evident that, initially, from the pilot symbols, one can 
obtain temporal estimated channel parameters corresponding to the pilot symbol positions. 
Then the channel parameters corresponding to the other positions are obtained by 
interpolation using a 2-D IFFT, filtering and a 2-D FFT. In the present estimator, the 2-D 
transforms and the multiplication are preferably performed using DSP chips, dedicated 
hardware, or the like. For instance, a single DSP chip can perform both 2-D FFTs. 

With the estimated channel parameters H m [n,k],a minimum mean square diversity 

combiner 128 can be used to combine the received signals x m [n,k] to estimate the combined 
signal as: 



The operation of the receiver 120 may now be explained with reference to the above 

equations and the channel estimator 140. As discussed above, the transmitter 100 and 

receiver 120 are synchronized so as to permit coherent signal detection. Therefore, the 

receiver knows the positions of the pilot symbols in any given OFDM block. During the first 

transmission of a sequence of OFDM blocks, the receiver uses only the known pilot symbol 

positions in the received signals x m [n,k], along with the known amplitudes and phases of the 

pilot symbols. During this first transmission, the receiver calculates a noisy estimate of the 
channel parameters in accordance with Eq. 2 to create H m [n,k] for just the known pilot 

symbol positions in the first sequence of OFDM blocks. A two-dimensional matrix 
G,[«, s « 2 ] comprising the calculated noisy channel parameters for the known pilot symbol 
positions, and zero-padded values for the remaining positions is formed. Matrix G,[w 1? « 2 ] is 
subjected to an IFFT, a diamond-shaped filter and an FFT to form the filtered channel 
estimation matrix G,[« is h 2 ] which comprises the initial noise-reduced channel estimates 



p 




(Eq. 3) 
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H m [n,k] received at the m antenna. The noise-reduced initial channel estimates are then 

used in Eq. 3 to estimate the combined signal yi[n,k] and the demodulated signal a y [n,k] . 

For a subsequent sequence of OFDM blocks, the receiver calculates a noisy enhanced 
channel estimate in accordance with Eq. 2 to create H m [n,k] using the m received signals 

5 x m [n 5 k] and the values of the estimated demodulated signal, rather than for only the n=s and 

k=s corresponding to the known pilot symbols. The same 2-D matrix is formed and 

subjected to an IFFT, a diamond filter and an inverse FFT to provide the noise-reduced 
enhanced channel estimates H m [n,k] for submission to Eq. 3, as before. 

After the initial channel estimate has been calculated, the positions of the pilot 

10 symbols do not solely determine the data that is used in calculating the enhanced channel 

estimates. Instead, the entire estimated demodulated signal is used in this calculation, and the 
channel estimate is optimum for that transmission. If, however, a transmission of a first set 
of sequences of OFDM blocks has concluded, and there is a time delay before a second set is 
transmitted, the previously estimated demodulated signal may not be suitable for use in the 

1 5 calculation of Eq. 2. In such case, for the second set of transmissions, it may be desirable to 
calculate a noise-reduced initial channel estimate using only the positions of the pilot 
symbols as explained above, and thereafter use the estimated demodulated signal for the 
subsequent noise-reduced enhanced channel estimates. 

In the present invention, at the OFDM transmitter, the pilot symbols are selectively 

20 inserted at predetermined tones in members of a sequence, or packet, of OFDM blocks. More 
particularly, the pilot symbol inserter 1 10 is configured to insert pilot symbols into successive 
OFDM blocks such that a predetermined pattern is produced on a two-dimensional grid 
formed from several rows of OFDM blocks, each row comprising the K tones of a 
corresponding block, each column comprising the signal at a particular tone for each of N 

25 OFDM blocks, with adjacent rows corresponding to adjacent OFDM blocks in the time 
domain. The placement of the pilot symbols is considered next. 

As discussed in P. Hoeher, S. Kaiser and P. Robertson, Two-dimensional pilot- 
symbol-aided channel estimation by Wiener filtering, Proc. of 1997 IEEE International 
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Conf. of Acoustics, Speech and Signal Processing, Munich, Germany, April 1997, pp 1845- 
1848 [1], and in P. Hoeher, S. Kaiser and P. Robertson, Pilot-symbol-aided channel 
estimation in time and frequency, Proc. of 1997 IEEE Global Telecomm. Conf: Comm. The. 
Mini-Conf, Phoenix, AZ, Nov. 1997, pp 90-96 [2], the grid density of the pilot symbols must 
5 satisfy the 2-D sampling theorem in order to recover channel parameters. For an OFDM 
signal, then, this means that the following must be satisfied: 

fa max T f M t < l/2;r max A/M y . < 1/2 (Eq 4) 

where fd max and T max are the maximum expected Doppler shift and maximum expected time 

10 delay, respectively; T f and Af represent the block length and the tone spacing, respectively, 
and M t and Mf denote the spaces between the time and frequency domains, respectively. It is 
suggested in the aforementioned articles that a twice oversampling rate for the pilot symbols 
be used to suppress channel noise and improve the estimation performance at the expense of 
reduced efficiency of the data transmission. 

15 In the aforementioned articles by P. Hoeher, S. Kaiser and P. Robertson, the authors 

studied parameter interpolation using Wiener filtering or the MMSE interpolation. These 
approaches to interpolation require information about the channel statistics. This contrasts 
with the present invention which provides for robust interpolation that is insensitive to 
channel statistics. The 2-D FFT and IFFT can be used to implement robust interpolation and 

20 reduce the computational complexity. However, in order to reduce edge effects, square-root 
raised cosine extensions are used on the edges of the 2-D signal, as disclosed in U.S. Patent 
No. 5,357,502, whose contents are also incorporated by reference herein. 

For wireless systems with packet data transmission, each packet may consist of only a 
few OFDM blocks. Therefore, the channel time-domain spectrum will be spread. For the 

25 same reason, the frequency-domain spectrum will also be spread due to a finite number of 
tones in each OFDM block. One may therefore increase the number of pilot symbols to 
reduce spectral aliasing at the expense of reduced efficiency of the data transmission. 
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However, for a given rate of pilot symbols, we may also improve performance by the choice 
of the pilot symbol grid shape. 

A rectangular grid shape in which pilot symbols are inserted periodically at the same 
tones in each of a number OFDM blocks is less preferable because such pattern does not 
5 efficiently capture variation in both the time delay and Doppler. In the present invention, one 
preferably uses a non-rectangular grid shape to take advantage of the two-dimensional 
filtering that is used in the estimator. The 2-D spectrum of a regular, non-rectangular pilot 
grid is discussed in the article R.M. Mersereau and T.C. Speake, The Processing of 
Periodically Sampled Multidimensional Signals, IEEE Trans, on Signal Processing, vol. 31, 

10 pp. 188-194, Feb, 1983 [3]. 

Fig. 4a shows a sample non-rectangular pilot symbol grid and Fig. 4b shows its 
corresponding 2-D spectrum contour. As seen in Fig. 4a, the pilot symbols are not 
necessarily inserted into adjacent rows, i.e., successive OFDM blocks. In addition, the pilot 
symbols in one OFDM block do not occupy the same tones as the pilot symbols in a nearby 

15 OFDM block, but rather are staggered relative thereto, thereby providing the non-rectangular 
pilot-symbol grid pattern. One skilled in the art will recognize from Fig. 4b that if a 
rectangular pilot symbol grid is used, there will be aliasing of its 2-D spectrum. 

Theoretical Performance 

20 As indicated above, with pilot symbols, temporal estimation of the corresponding 

parameters can be obtained. The channel parameters corresponding to the regular symbols 
(unknown data) can be obtained by interpolation. The optimum interpolation is determined 
by the channel statistics that are usually difficult to obtain since they depend on the 
environment and vehicle speed. We now demonstrate that the present invention results in a 

25 robust interpolator. 

MMSE Interpolation. The channel parameter H[p^p 2 ] , which is a two-dimensional 
sequence, will be written as . H[p] for simplicity, where pk(p x >p 2 ) T denotes a 
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2-element integer vector. From now on, the subscript m for antenna is omitted since the 
channel estimators are the same for different antennas. From the pilot symbols with any 
regular grid, we can get temporal estimation of the channel parameters, 



H[Vp) = H[Vp] + <D[Vp), (Eq.5) 
for all p, where V is a non-singular 2x2 matrix that is determined by the pilot symbol grid. 
For the grid shown in Figure 3(a), 



V = 



2 5 
2 0 



(Eq. 6) 



1 0 For any 2-element integer vector n = (ni.nif, there is a unique pair of p and r / (where / is a 
set of |det(F)| 2-element integer vectors determined by V), such that 

n = Vp + r (Eq- 7) 



From the temporal estimation, H[n], for any n, can be estimated by 

15 

H[n] = ^c[r + V(p-p 0 )]H[Vp 0 l (Eq. 8) 

Pa 

where the c[n]'s are two-dimensional coefficients. Using the orthogonality principle, the 
coefficients for the MMSE estimator are determined by solving 

20 £{(^[n]-Z4r + K(p-pJMFpJ^[K Pl ]} = 0 5 (Eq. 9) 

for all pi, or 

r[Vp + r] = ^c[r + V(p-p o )MVp 0 ) + pS[p 0 ]), (Eq. 10) 
p. 
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where 

p = E{ai[Vp]\ 2 } 



Denote 

R(r,aj)^r[Vp + T]e- J ^, (Eq.ll) 



CM^Wp + rV^"", (Eq. 12) 

P 

where co = (0)1,0)2/. Then, (10) can be expressed in the frequency domain as 

10 R(r,a)) = C(r,a))(R(0,o)) + p). (Eq. 13) 

Therefore, 

R(0, co) + p 

Let R(co) be the 2-D Fourier transform of r[n] defined as 
15 *(fl» = £r[n>-** (Eq. 15) 

n 

Then, we have 

/?(r,^)-^r[Kp-hrK^ p 
p 

= I{^T £ f/(^* W) ^}^ (Eq. 16) 

1 



2 £ £^K){X^ ( "^ )r,,, }e y ' fr rf<^, 

p 



(2*) 
20 

with =(#",^2)- Since, 



Xe-^^ir) 2 ^^^-^-^). (Eq. 17) 

P P 

Then, 

13 



R(r, a) = —- - r r R{(o 0 ){2n) 2 2 S{V T a> - V T co„ - In v).e^ dco" x da>\ 



PeZ 1 



(Eq. 



with U=V X . 



From (18), we have that 



re/ 



= V_L_ y ^(^- 2 ^P)^"" ((/P)r)f e ^ (Eq 19) 



It can be proven from the property of /that 



e J 



Consequently, 



re/ [ 0, otherwise. 



CW = ^f- (Eq.21) 
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If the pilot symbols are dense enough to make R{r,co) alias-free, then R(0,co) = 
within the non-zero support of R(co) and 



|det(F)| 



|det(F)| +P 



which implies that the optimum parameter interpolator can be determined once the channel 
statistics are known. 

Mean-Square Error (MSE). The MSE for the optimum interpolation can be found by 

mse 0 =E{\h[Vp + r)- £ c[T,V(p- Po )}h[Vp 0 ]\ 2 } 

= E{h[Vp + r] - X c t r . K(p - p 0 )]h[Vp 0 ]h'[Vp + r]} 

(Eq. 23) 

= r[0]- Y c l r > V (Po)V[Vp a +r] 

P„6Z 2 

= 1 — L_ r rJ^c^L^ 

If the pilot symbols satisfy the alias-free condition, then |/?(r,(y)| 2 = |i?(0,«)| 2 . 
Furthermore, there are |det(F)| replicas of rj - ^^ within [-;r,7r]x[-;r,7r]. Let the 

non-zero support of R(m) be D with area S, then 



|det(FQ| ff \R(0,a>)\ 2 j 

mse„ =1- J r 1 J 1 — da>,de)i 

(2tt) 2 ii ^R{0,w) + p 1 " 



R(a>) 



(2;r) 2 J * 1 2 (2;r) 2 J -b AQa) 1 2 

|det(F)| + P 



(2^) 2 J ^ *(*) , 



|det(0| 
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If R(ta) is ideal low-pass, that is 



(2/r) 2 
0, 



otherwise, 



then 



5|det(F)| 



+ J0 



(Eq. 25) 



(Eq. 26) 



For any ^?f<y) with non-zero support D and satisfying 



10 we have 



{2n)' 



■\^R{ca)da x da> 1 =\ (Eq. 27) 



S 1 f f K(0)) 



R(a>) 



-dco.da>. 



|det(F)| 
^r\] D R{o>)d(O x do} 2 



(2^ 1 rr_^L 

S J *>|de«7)| 1 2 ^ 



(Eq. 28) 



(2*) z 
S|det(7)| 



where we have used the fact that is a convex function for x > 0. The inequality (28) 

x + p 

implies that of all R(co) with non-zero support D, the channel with the ideal bandlimited 2-D 
spectrum has the worst performance. 
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Robust interpolation. The MMSE interpolation derived in the previous section 
requires channel statistics that are usually unknown since it changes with the environment 
and vehicle speed. Hence, a robust interpolator that is independent of the channel statistics is 
desired. 

For an interpolator matching the 2-D spectrum R(co) and a channel with 2-D 
spectrum R{a>) , the MSE of the interpolator will be 



<r 4 

mse = 



a" rc R(o))-R(co) . , _ _ 

|det(F)| + P) 



10 

The first term in the above equation is the MSE increment due to the mismatch of the channel 
statistics and the estimator. However, if R{co) = R r {co) then mse = mse 0 and the performance 

of the interpolator does not depend on the channel statistics. We call the interpolator that 
matches R r {o)) a robust interpolator. The robust interpolator can be shown as in Figure 3. 

1 5 Since the 2-D spectrum is usually of a diamond shape, a diamond shape filter is used to filter 
out the additive noise. 

Enhanced parameter estimation. If an M times oversampling rate is used for the pilot 
symbols, then S = {In) 1 IM and from Eq. (26), the MSE of the robust channel estimator will 
20 be 

w ^ = Trb*T7 (Eq 3 °) 

M + p M 

An enhanced channel parameter estimation algorithm as discussed in Y. (Geoffrey) Li 
an N. Sollenberger, Interference suppression in OFDM systems using adaptive antenna 
25 arrays, 1998 Global Telecomm. Conf: Comm. The. Mini-Conf, Sydney, Australia, Nov. 

1998 [4], can be used here to improve the performance of the channel parameter estimation. 
Based on the initial parameter estimation using the pilot symbols, the transmitted symbols 
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can be coherently detected, and therefore the temporal estimation for all Ids and j/s can be 
obtained with the help of the detected symbols. Better parameter estimation can then be 
obtained from the temporal parameter estimator. If all detected symbols are correct, then 
from (26) the MSE of the enhanced parameter estimator is 

mse = ^ » P - . (Eq. 31) 

6 (2tt) 2 (2tt) 2 14 ' 

- — — + p - — — 



Since — is usually much larger than M, the enhanced estimation can significantly 

S 

improve the system performance even though there are symbol detection errors and edge 
effects, as demonstrated by the simulations described below. 

10 The simulation uses channels having one of three delay profiles: (1) the two-ray delay 

profile of Y. (Geoffrey) Li, J.H Winters and N.R. Sollenberger, Spatial-temporal equalization 
for IS- 136 TDMA systems with rapid dispersive fading and co-channel interference, IEEE 
Trans, on Veh. Tech., to be published [5]; (2) typical-urban (TU) delay profile and (3) hilly- 
terrain (HT) delay profile, both discussed in Y. (Geoffrey) Li, N. Seshadri and S. 

1 5 Ariyavisitakul, Transmitter diversity in OFDM systems with dispersive fading channels, 

presented at 1998 IEEE Global Telecomm. Conf., Sydney Australia, November 1998 [6] and 
also in Y. (Geoffrey) Li, N. Seshadri and S. Ariyavisitakul, Channel estimation for 
transmitter diversity in OFDM systems with mobile wireless channels, IEEE J-SAC Wireless 
Comm. Series, to be published [7]. Each of the three delay profiles was considered with both 

20 40 Hz and 200 Hz Doppler frequencies. The channels corresponding to different receivers 
have the same statistics and two antennas are used for receiver diversity. 

The parameters for the OFDM system are the same as those in Y. (Geoffrey) Li, L. 
Cimini and N. Sollenberger, Robust Channel Estimation for OFDM systems with rapid 
dispersive fading channels, IEEE Trans, on Comm., vol. 46, pp. 902-915, July 1998 [8]. To 

25 construct an OFDM signal, assume the entire channel bandwidth, 800 kHz, is divided into 
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128 tones. The 4 tones on each end are used as guard tones, and the rest (120 tones) are used 
to transmit data. To make the tones orthogonal to each other, the symbol duration is 160 josec. 
An additional 40 jasec guard interval is used to provide protection from intersymbol 
interference due to channel multipath delay spread. This results in a total block length T = 
5 200 (isec and a subchannel symbol rate rb = kbaud. 

PSK modulation with coherent demodulation is used. As in reference [8], a (40,20) 
R-S code, with each code symbol consisting of 3 QPSK symbols grouped in frequency, is 
used in the system. The R-S decoder erases 10 symbols based on the signal strength and 
corrects 5 additional random errors. The system is in packet transmission mode with each 
10 packet containing 10 OFDM blocks. Pilot symbols are inserted in the system at a 10% rate 
with the grid shown in Fig. 4a. Hence, the simulated system can transmit data at 1.08 
Mbits/sec before decoding, or 540 kbits/sec after decoding, over an 800 kHz channel. To get 
insight into the average behavior of the channel estimator, the performance was averaged 
over 10,000 OFDM blocks. 

15 

Effects of different references on system performance 

Figs. 5a and 5b show the effects of different references on the word error rate (WER) 
and the mean squared error (MSE) of the OFDM system with the TU channel and 40 Hz 
Doppler frequency. From these figures, when 10% pilot symbols are used for initial 

20 estimation, the MSE of the estimated parameters is almost the same as the system's SNR. 

The required SNR's for 10% and 1% WER's are about 8 dB and 1 1 .2 dB, respectively. When 
the decisioned (sliced) symbols along with the pilot symbols are used to enhance the channel 
parameter estimation, the performance of parameter estimation can be improved significantly, 
and the MSE of the enhanced estimator is reduced from -10 dB to -15.5 dB when SNR=10 

25 dB. The required SNR's for 10% and 1% WER's are improved by about 1.4 dB. If the 

undecoded/decoded dual mode reference introduced in reference [8] is used, another 0.3 dB 
SNR improvement can be obtained. However, when an ideal reference is used for channel 
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parameter estimation, the system performance is much better than that of initial or enhanced 
channel parameter estimation. 

Parameter estimation performance under different environments 
5 Figs. 6a, 6b and 7a, 7b show the performance of the initial and enhanced parameter 

estimators for channels with the two-ray (having a delay span of 5 |osec), TU, and HT delay 
profiles and 40 Hz and 200 Hz Doppler frequencies, respectively. From the figures, the 
system with a pilot-symbol-aided estimator is very robust to Doppler shift, and there are 
almost no performance degradations for the two-ray and TU channels and only 0.3 dB SNR 
10 degradation for the HT channel if fa is increased from 40 Hz to 200 Hz. The OFDM system 
for the two-ray and TU channels has almost the same performance, which is better than the 
system for the HT channel. The required SNR for a 10% WER is 7 dB for the two-ray and 
TU channels and about 8 dB for the HT channel. 

1 5 Comparison with decision-directed estimator under noise-limited performance 

Figs. 8a and 8b compare the pilot-symbol-aided and the decision-directed (see 
Reference [8]) estimators for an OFDM system with the TU and HT channels and 40 Hz and 
200 Hz Doppler frequencies. From Fig. 8a, for the TU channel with fd = 200 Hz, the pilot- 
symbol-aided and the decision-directed estimators have the same performance and the 

20 required SNR's for 10% and 1% WER's are about 7 dB and 10 dB, respectively. However, 
for the TU channel with f d = 200 Hz, the performance of the decision-directed estimator is 
degraded significantly, and the required SNR for a 10% WER is as large as 9.3 dB while it is 
only about 7.3 dB for the pilot-symbol-aided estimator. Fig. 8b demonstrates the WER's of 
system for the HT channel. It shows that the decision-directed estimator is 1 dB better than 

25 the pilot-symbol-aided estimator for lower Doppler frequency (f d = 40 Hz); however, the 
pilot-symbol-aided estimator is much better for higher Doppler frequency (fa = 200 Hz). 

From the foregoing simulation for evaluating pilot-symbol-aided channel parameter 
estimation for OFDM systems, it appears that robust interpolation using the 2-D FFT and 

20 



IFFT, and enhanced parameter estimation can significantly improve the performance of 
OFDM in noise-limited wireless environments. Compared with decision-directed channel 
parameter estimator, the pilot-symbol-aided estimator is highly robust to Doppler, even 
though it has a slight performance degradation for systems with lower Doppler frequencies. 
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